home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / clahqr.z / clahqr
Encoding:
Text File  |  2002-10-03  |  5.0 KB  |  133 lines
  1.  
  2.  
  3.  
  4. CCCCLLLLAAAAHHHHQQQQRRRR((((3333SSSS))))                                                          CCCCLLLLAAAAHHHHQQQQRRRR((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      CLAHQR - i an auxiliary routine called by CHSEQR to update the
  10.      eigenvalues and Schur decomposition already computed by CHSEQR, by
  11.      dealing with the Hessenberg submatrix in rows and columns ILO to IHI
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE CLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, W, ILOZ, IHIZ, Z,
  15.                         LDZ, INFO )
  16.  
  17.          LOGICAL        WANTT, WANTZ
  18.  
  19.          INTEGER        IHI, IHIZ, ILO, ILOZ, INFO, LDH, LDZ, N
  20.  
  21.          COMPLEX        H( LDH, * ), W( * ), Z( LDZ, * )
  22.  
  23. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  24.      These routines are part of the SCSL Scientific Library and can be loaded
  25.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  26.      directs the linker to use the multi-processor version of the library.
  27.  
  28.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  29.      4 bytes (32 bits). Another version of SCSL is available in which integers
  30.      are 8 bytes (64 bits).  This version allows the user access to larger
  31.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  32.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  33.      only one of the two versions; 4-byte integer and 8-byte integer library
  34.      calls cannot be mixed.
  35.  
  36. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  37.      CLAHQR is an auxiliary routine called by CHSEQR to update the eigenvalues
  38.      and Schur decomposition already computed by CHSEQR, by dealing with the
  39.      Hessenberg submatrix in rows and columns ILO to IHI.
  40.  
  41. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  42.      WANTT   (input) LOGICAL
  43.              = .TRUE. : the full Schur form T is required;
  44.              = .FALSE.: only eigenvalues are required.
  45.  
  46.      WANTZ   (input) LOGICAL
  47.              = .TRUE. : the matrix of Schur vectors Z is required;
  48.              = .FALSE.: Schur vectors are not required.
  49.  
  50.      N       (input) INTEGER
  51.              The order of the matrix H.  N >= 0.
  52.  
  53.      ILO     (input) INTEGER
  54.              IHI     (input) INTEGER It is assumed that H is already upper
  55.              triangular in rows and columns IHI+1:N, and that H(ILO,ILO-1) = 0
  56.              (unless ILO = 1).  CLAHQR works primarily with the Hessenberg
  57.              submatrix in rows and columns ILO to IHI, but applies
  58.              transformations to all of H if WANTT is .TRUE..  1 <= ILO <=
  59.              max(1,IHI); IHI <= N.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. CCCCLLLLAAAAHHHHQQQQRRRR((((3333SSSS))))                                                          CCCCLLLLAAAAHHHHQQQQRRRR((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      H       (input/output) COMPLEX array, dimension (LDH,N)
  75.              On entry, the upper Hessenberg matrix H.  On exit, if WANTT is
  76.              .TRUE., H is upper triangular in rows and columns ILO:IHI, with
  77.              any 2-by-2 diagonal blocks in standard form. If WANTT is .FALSE.,
  78.              the contents of H are unspecified on exit.
  79.  
  80.      LDH     (input) INTEGER
  81.              The leading dimension of the array H. LDH >= max(1,N).
  82.  
  83.      W       (output) COMPLEX array, dimension (N)
  84.              The computed eigenvalues ILO to IHI are stored in the
  85.              corresponding elements of W. If WANTT is .TRUE., the eigenvalues
  86.              are stored in the same order as on the diagonal of the Schur form
  87.              returned in H, with W(i) = H(i,i).
  88.  
  89.      ILOZ    (input) INTEGER
  90.              IHIZ    (input) INTEGER Specify the rows of Z to which
  91.              transformations must be applied if WANTZ is .TRUE..  1 <= ILOZ <=
  92.              ILO; IHI <= IHIZ <= N.
  93.  
  94.      Z       (input/output) COMPLEX array, dimension (LDZ,N)
  95.              If WANTZ is .TRUE., on entry Z must contain the current matrix Z
  96.              of transformations accumulated by CHSEQR, and on exit Z has been
  97.              updated; transformations are applied only to the submatrix
  98.              Z(ILOZ:IHIZ,ILO:IHI).  If WANTZ is .FALSE., Z is not referenced.
  99.  
  100.      LDZ     (input) INTEGER
  101.              The leading dimension of the array Z. LDZ >= max(1,N).
  102.  
  103.      INFO    (output) INTEGER
  104.              = 0: successful exit
  105.              > 0: if INFO = i, CLAHQR failed to compute all the eigenvalues
  106.              ILO to IHI in a total of 30*(IHI-ILO+1) iterations; elements
  107.              i+1:ihi of W contain those eigenvalues which have been
  108.              successfully computed.
  109.  
  110. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  111.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  112.  
  113.      This man page is available only online.
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.